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Abstract 

We present exact calculations of the zero-temperature partition function (chromatic poly- 
nomial) P for the g-state Potts antiferromagnet on triangular lattice strips of arbitrarily great 
length Lx vertices and of width Ly = 3 vertices and, in the —>■ oo limit, the exponent of 
the ground-state entropy, W = e^°/^^ . The strips considered, with their boundary conditions 
{BC) are (a) {FBCy, PBC^) = cychc, (b) {FBCy, TPBC^) = M5bius, (c) {PBCy, PBC^) = 
toroidal, and (d) {PBCy, TPBC^) = Klein bottle, where F, P, and TP denote free, periodic, 
and twisted periodic. Exact calculations of P and W are also given for wider strips, includ- 
ing (e) cyclic, Ly = 4, and (f) {PBCy, FBCx) = cylindrical, Ly = 5, 6. Several interesting 
features are found, including the presence of terms in P proportional to cos(27rLj./3) for case 
(c). The continuous locus of points B where W is nonanalytic in the q plane is discussed 
for each case and a comparative discussion is given of the respective loci B for families with 
different boundary conditions. Numerical values of W are given for infinite-length strips of 
various widths and are shown to approach values for the 2D lattice rapidly. A remark is also 
made concerning a zero-free region for chromatic zeros. 
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1 Introduction 



The g-state Potts antiferromagnet (AF) [|l], |^ exhibits nonzero ground state entropy, 5*0 > 
(without frustration) for sufficiently large g on a given lattice A or, more generally, on a 
graph G. This is equivalent to a ground state degeneracy per site W > 1, since 5*0 = 
ks^^W . Such nonzero ground state entropy is important as an exception to the third law 
of thermodynamics. There is a close connection with graph theory here, since the zero- 
temperature partition function of the above-mentioned g-state Potts antiferromagnet on a 
graph G satisfies 

Z{G,q,T = Q)pAF = P{G,q) (1.1) 

where P(G, q) is the chromatic polynomial expressing the number of ways of coloring the 
vertices of the graph G with q colors such that no two adjacent vertices have the same color 
(for reviews, see 0-[§l). The minimum number of colors necessary for such a coloring of G 
is called the chromatic number, xiG). ThusQ 

W{{G},q)= hm P(G,g)i/" (1.2) 

71— >00 

where n = v{G) is the number of vertices of G and {G} = limn^^G. Since P{G,q) is 
a polynomial, one can generalize q from Z+ to C. The zeros of P{G, q) in the complex 
q plane are called chromatic zeros; a subset of these may form an accumulation set in 
the n — s> oo limit, denoted B, which is the continuous locus of points where W{{G},q) is 
nonanalytic. The maximal region in the complex q plane to which one can analytically 
continue the function W{{G},q) from physical values where there is nonzero ground state 
entropy is denoted Ri. The maximal value of q where B intersects the (positive) real axis is 
labelled qc{{G}). This point is important since it separates the interval q > qc{{G}) on the 
positive real q axis where the Potts model (with q extended from Z+ to M) exhibits nonzero 
ground state entropy (which increases with q, asymptotically approaching 5*0 = feeing for 
large q, and which for a regular lattice A can be calculated approximately via large-g series 
expansions) from the interval < g < qc{{G}) in which 5*0 has a different analytic form. 

^At certain special points Qs (typically — 0, 1, .., x(G)), one has the noncommutativity of limits 
lim^^q^ lini„^oo P(G, g)^/" ^ Hm„^oo limq_^^ P(G', g)^/", and hence it is necessary to specify the order 
of the limits in the definition of W^dG}, Qs) 0. We use the first order of limits here; this has the advantage 
of removing certain isolated discontinuities in W. 

^Although it does not happen in the cases considered here, for some families of graphs B may be null, 
and W may also be nonanalytic at certain discrete points. 
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Figure 1: Illustrative strip graphs of the triangular lattice with width Ly = 3 and length ^ 4: hav- 
ing the following boundary conditions: (a) {FBCy^PBC^) = cyclic, (b) {FBCy,TPBCx) = Mobius, (c) 
{PBCy,PBC^) = toroidal, and (d) {PBCy,TPBC^) = Klein bottle. 

In the present work we report exact calculations of P{G, q), W{{G}, g), and B for families 
of triangular lattice strips of fixed width Ly — 3 vertices and arbitrary length — m vertices. 
The longitudinal and transverse directions are taken to be x and y). In Fig. 1 we display 
some illustrative examples. A generic form for chromatic polynomials for recursively defined 
families of graphs, of which strip graphs Gg are special cases, is 

P((aUg) = ^c,(g)(A,(g)r (1-3) 

where Cj{q) and the Nx terms Xj{q) depend on the type of strip graph Gg but are independent 
of m. 
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Let us comment on the motivations for the current study and how it extends previous 
work. Calculations of the chromatic polynomials for the cyclic and Mobius strip graphs 
of the square lattice were carried out for Lj^ = 2 in [§] (see also |10|)- Subsequently, 



these calculations were extended to the width Ly = 3 for the cyclic ||TT| and Mobius |jT2| 
square lattice strips. After studies of chromatic zeros in [jl3|, [T^ and, for Ly = 2 strips in 



T5| , [T^, the W and B defined in the infinite-length limit were determined for Ly 



in 0] and for = 3 in |]TT], |T^. Similar calculations were carried out for strips of various 
lattices with widths up to L^, = 4 and free boundary conditions in [|l^, |18|, |19|, |2^, and 
for cyclic and Mobius graphs involving homeomorphic expansions of square strips and 
cyclic polygon chains ^2,\. An important question concerns the effect of boundary conditions 
(BC's), and hence graph topology, on P, W, and B. We use the symbols FBCy and PBCy 
for free and periodic transverse boundary conditions and FBCx, PBC^^, and TPBC^; for free, 
periodic, and twisted periodic longitudinal boundary conditions. The term "twisted" means 
that the longitudinal ends of the strip are identified with reversed orientation. These strip 
graphs can be embedded on surfaces with the following topologies: (i) (FBCy,FBCa;): open 
strip; (ii) (PBCj,,FBC2;): cylindrical; (iii) (FBCj/,PBCa;): cylindrical, denoted cychc here; 
(iv) (FBCj,,TPBa): M5bius; (v) (PBCy,PBa.): torus; and (vi) (PBCj^TPBa): Klein 
bottle (We recall that that unlike graphs of type (i)-(v), the Klein bottle surface cannot 
be embedded without self-intersection in M^.) 

Several specific effects of different boundary conditions have been found (for a general 
discussion, see p3| , p5|). First, from comparisons of exact calculations of W for the infinite- 
length limit of lattice strips with boundary conditions of types (i)-(iv), it was observed 
that those with a global circuit^ lead to a locus B that separates the q plane into different 
regions, and it was inferred that the presence of global circuits is a sufficient condition for 
B to have this property 0, 0, ^ |ll], |12|-Q For strip graphs of regular lattices, this is 
equivalent to having PBCx or TPBC^, i.e., periodic or twisted periodic boundary conditions 
in the direction in which the strip length goes to infinity as — > oo. Further support for 
this inference was adduced from the calculation of P, W, and B for the Ly = 3 strip of 
the square lattice with torus and Klein bottle boundary conditions [^]. Applications of a 
coloring matrix method [|l^ for these calculations has been discussed in |p7| , |28| . For the 
(homogeneous) strip graphs with (FBCy,FBC2:) in [l^, |l9l, B consists of arcs that do not 

■^These BC's can all be implemented in a manner that is uniform in the length L^; the case (vii) 
(TPBCj,,TPBCa;) with the topology of the projective plane requires different identifications as varies 
and will not be considered here. 

"^A global circuit is a route following a lattice direction which has the topology of and a length ig,c. 
that goes to infinity as n — s- oo. 

^This is not a necessary condition, as was shown in 
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separate regions of the q plane. As the width of the strip increases, these arcs tend to 
elongate, and their ends tend to move toward each other, thereby suggesting that if one 
considered the sequence of strip graphs of this type with width Ly (having taken the limit 
Lj. — i> oo to obtain a locus B for each member of this sequence), then in the limit Ly oo, 
the arcs would close to form a B that separated the q plane into two or more regions and 
passed through g = 0. The interesting feature of the families of graphs with global circuits is 
that when the graphs contain a global circuit, this separation of the q plane into regions (with 
B passing through g = 0) already occurs for finite Ly. This means that the W functions of 
the infinite-length limit of these graphs with cyclic and twisted cyclic longitudinal boundary 
conditions already exhibit a feature which is expected to occur for the B for the W function 
of the full two-dimensional lattice. This expectation is supported by the calculation of W and 



B for the 2D triangular lattice with cylindrical boundary conditions by Baxter ||2^. Thus, 
although calculations of chromatic polynomials for lattice strips with global circuits are, in 
general, more difficult than for the corresponding strips with free boundary conditions, such 
calculations are worthwhile, since in the L^ ^ oo limit, the resultant locus B embodies the 
analytic property of B expected for the full 2D lattice, viz., that it crosses the real axis at 
least at g = and a maximal point qc, and separates the q plane into different regions. 

Second, it has been found that for a given type of strip graph Gs with FBC^, the chromatic 
polynomials for PBC^. and TPBC^ boundary conditions (i.e., cyclic and Mobius strips) have 
the same Xj, although in general different Cj. It follows that the loci B are the same for these 
two different longitudinal boundary condition choices |2l|, 0, |12|. In the case of PBCy, 



the reversal of orientation involved in going from PBC^ to TPBC^. longitudinal boundary 
conditions (i.e. from torus to Klein bottle topology) can lead to the removal of some of the 
Aj's that were present; i.e., P for the (PBCy,TPBC2:) strip may involve only a subset of 
the Aj's that are present for the (PBCj,,PBC2:) strip. For example, for the Ly = 3 strips of 
the square lattice with (PBC^, PBC^;) boundary conditions, there are A^a = 8 A/s, but for 
the strip with (PBC^, TPBC^;) boundary conditions only a subset of iV^ = 5 of these terms 
occurs in P None of the three Aj's that are absent from P in the TPBC^; case is leading, 
so that B is the same for both of these families. 

Third, for a given type of strip graph Gs containing a global circuit, it has been found 
that, in the infinite- length limit where the locus B is defined, it not only passes through the 
origin, q = 0, but always crosses the positive real axis at one or more points, the maximal 
one being denoted qc{{G}), as mentioned above. In contrast, for strip graphs that do not 
contain global circuits, B may not cross the real axis. For example, for the (infinite-length 
limit of the) strip of the triangular lattice with free transverse and longitudinal boundary 
conditions, {FBCy, FBG^), the locus B crosses the real axis for width Ly = 3 but does not 
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cross it for Ly = 4 (see Fig. 5 of [|T3l)- For square strips with {FBCy, (T)PBCx) boundary 
conditions, it has been found that qc{{G}) is a nondecreasing function of Ly for the cases 



studied so far, namely Qc = 2 for Ly = 1 and Ly = 2, and Qc — 2.33654 for Ly = 3 [|lT|, 0, P3| . 
For the strips of the triangular lattice with {FBCy, {T)PBCx) boundary conditions, we only 
have one width for which B has been determined, namely, Ly = 2, for which Qc = 3. It is thus 
of interest to explore what the value of qc is for cyclic and Mobius strips of the triangular 
lattice with larger widths. For the Ly = 3 torus and Klein bottle strips of the square lattice, 
very interestingly, = 3, which is the same value as for the 2D square lattice ||2^. Clearly 
it is worthwhile to investigate what the value of qc is for torus and Klein bottle strips of the 
triangular lattice, and we shall do this here. 

Fourth, for the strip graph {Gs)m with a given type of transverse boundary conditions 
BCj^, the chromatic poljTiomial for PBCa; has a larger value of than the chromatic poly- 
nomial for FBCx, and the corresponding loci B are different. 

Fifth, one may ask how the W functions in region Ri compare for different boundary 



conditions. It has been found [^, |22|, |TT], |T^, ^ that, for a given type of strip graph 
Gs, in the region Ri (with its left-hand boundary on the positive real q axis at qc) being 
defined for the PBCa; boundary conditions, the W function is the same for FBCa;, PBC^^, 
and TPBC^. This includes the region of real q greater than this value of qc, so that this 
result is somewhat reminiscent of the statement of the existence of the thermodynamic limit 
in statistical mechanics, i.e. the fact that in the disordered phase of a statistical mechanical 
system the thermodynamic functions are independent of the boundary conditions and if an 
ordered phase exists, then the role of the boundary conditions is only to set a preferred 
direction for the symmetry-breaking order parameter. 

Our new results on P, W, and B for the Ly = 3 strip of the triangular lattice with the 
various boundary conditions of type (iii)-(vi) give insight concerning these five items, as do 
our additional results for Ly = 4 cyclic strips and Ly = 5,6 cylindrical strips. Calculations 



of these quantities for the Ly = 2 cyclic and Mobius triangular strips were presented in [|lT 



An early study of P for the Ly = 2 strip is in and a coloring matrix method to obtain 
P has been given in |2^. P, W, and B have been calculated for triangular lattice strips with 
free boundary conditions (i) for Ly up to 4 [|l^. (Some related calculations including work 
on cases with noncompact B, are |jl6|, 0], |3ll 



For strips with global circuits, and L^ above the lowest values involving degenerate cases, 
the triangular lattice strips with boundary conditions of type (i)-(vi) have n = LyL^ vertices; 

^We note the following corrections in Sand's result for P for the Ly = 2 cyclic strip of the triangular 
lattice, given on p. 88 of the coefficients of the last two terms should be reversed in sign to read {z = q 



in his notation) (z - I) [(1/2) (5 - 2z + ^9 - 4z)]™ + [(1/2) (5 - 2z - y/9 - 4z)]" , as in eq. (5.23) of p 
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the cyclic and Mobius triangular strips have e = {3Ly — 2)Lx edges while the torus and Klein 
bottle triangular strips have e = SL^L^. edges. Define Nt to be the number of triangles 
in the strips (i.e., excluding the triangles forming the cross sections of the strips in the 
case Ly = 3); then, again for and Ly large enough to avoid degenerate cases, the cyclic 
and Mobius strips have Nt = 2{Ly — l)Lx such triangles, while the torus and Klein bottle 
strips have 2LyLx such triangles. For strips of type (i), e = SL^Ly — 2Lx — 2Ly + 1 and 
Nt = 2{Ly — l){Lx — 1) while for strips of type (ii), e = SL^Ly — 2Ly and Nt = 2Ly{Lx — 1). 

Again, for = m and Ly large enough to avoid degenerate cases, the chromatic number 
for the triangular lattice strips with cyclic or torus boundary conditions is 



x{trz{Ly)m, FBCy, PBC,) = x{tri{Ly)„,, PBCy, PBC,) 



3 if m = mod 3 

4 if m = 1 or 2 mod 3 

(1.4) 

and for the triangular lattice strips with Mobius or Klein bottle boundary conditions. 



x{tri{LyU, FBCy, TPBC,) = x{tri{Ly)„., PBCy, TPBC,) = 4 V m > 3 . (1.5) 
For strips with (FBCy, FBC^) or {PBCy, FBC^), x = 3. 

2 Ly = 3 Cyclic and Mobius Triangular Strips 

We calculate the chromatic polynomials by iterated use of the deletion-contraction theorem, 
via a generating function approach [l^, |19|. In general, for a strip graph of type Gg we have 



r{Gs,q,x)= P{{Gs)m,q)x"' (2.1) 



m=mo 



where, as before [|TT|, we take mo = 2 for strips with periodic or twisted periodic longitudinal 
boundary conditions. The generating functions T{Gs, q, x) are rational functions of the form 

r(G..,.x) = ^J^ (2,2) 

with 

M{Gs,q,x) = Y,AGsM)^' (2-3) 

3=0 

and 

V[G„q,x) = l + Y,bG.My (2.4) 
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where the Ac^^i and bc^^i are polynomials in q (with no common factors) and 



dj^ = deg^{J\f) 



(2.5) 



and 

dv = deg^iV) . 
The denominator V can be written in factorized form as 



A 



Gs,j 



X] 



(2.6) 



(2.7) 



V{Gs,q,x) 

Equivalently, the XgsJ ^^re roots of the equation 

^"-ViGs, g, I/O = + X] &G.,,e'"-^' • (2.8) 
From the generating function, we calculate the chromatic polynomials in the form of eq. 



|1.3D by the use of the general formula, eq. (2.14), in [19 



d-D r d^r 

j=i Ls=o 



n 



\ rn—rno 



(2.9) 



P{Gm,q) 

(where the conventions in eq. (|2.1|) were chosen such that mg = in [|19[)- Note that some 
of the coefficients cq^.j may vanish, so that not all of the Ag^j's in V{Gs, q, x) contribute to 
the sum in (|1.3| ) ||11|| . The formula (|1.3|) for the chromatic polynomial has the virtue of being 
a closed-form expression. However, for a given type of strip graph, for values of the width 
Ly greater than the first one or two values, the set of Aj's include nonpolynomial algebraic 
roots. These occur, of course, as symmetric functions of the roots, so that, by a theorem 
on symmetric functions of roots of algebraic equations, their resultant contributions to P 
are expressible in terms of the coefficients of these equations (which are polynomials in q) 



and hence are polynomials in q [0]. However, in calculating chromatic polynomials and the 
resultant chromatic zeros for large values of = m, it can be convenient to use directly 
the expression for the generating function, since both the numerator and denominator of 
this function, eqs. (|2.3|), ( |2.4| ), are polynomials in q so that one does not have to rely 
upon cancellations of nonpolynomial algebraic expressions at intermediate stages in the 
calculation. 

The coefficients Ctri{Ly),j that enter into the expressions for the chromatic polynomial 
( |1.3|) for the cyclic triangular (t) strip of width Ly are (like those for the cyclic square strip) 
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certain polynomials that we denote c^'^\ given by |^ 



n(^- 



c^''^ = I l[q- qd.k) 
k=l 



where 



grffc = 2 + 2cos( — ; ) for k = l,2,...d 

^' \2d+lJ 



(2.10) 



(2.11) 



with < d < Ly. We list below the specific c^'^^'s that appear in our results for the Ly = 3 
and Ly = 4 cyclic strips of the triangular lattice 



c(°) = 1 , c^'> =q-l , d''> =q'-3q + l , 



(2) - ^2 



-,(3) = ^3 _ c- 2 



- + 6g - 1 



(2.12) 
(2.13) 



and 



c(^) = (g-l)(g3_gg2^g^_^^) _ ^^2.14) 

Where the coefficient c''°^ appears in chromatic polynomials, we shall simply write it as unity. 



For the cyclic L„ = 3 triangular strip of length = m, we obtain the result 
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P(trz(3 X m),FBCy,PBC,,q) = J] Ct3j(At3j) 

c(3)(_l)- + c(2) \{q - 2r + (A«,3)™ + (Ai3,4)' 



+ c 



(1) 



-(g-2)(g-3)]™ + (Ai3,6)™ + (A«,7)" + (A«,8)™ + (A^.g)'" + (Ais.io) 



(2.15) 



where the short notation t3 is used in subscripts for tri{Ly = 3). The labelling of and 
Xt3j in eq. ( P.15D is consecutive, so that Ct3,i = c^^\ Xt3,i = -1; Ct3,2 = c^^^, Af3,2 = g - 2, 
and so forth: 



A 



t3,(3,4) 



^(-7 + 2g± v/25-8g 



At3 7, j = 6, 7, 8 are the roots of the cubic equation 



where 



+ bt3,i^^ + bt3,2^ + bt3,3 = 

bt3,i = 2g' - 12g + 19 



(2.16) 

(2.17) 
(2.18) 
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6t3,2 = (g - 2)(g3 - 9g2 + 28q - 29) (2.19) 
bt3,3 = -{q-2)\q-3) , (2.20) 

and 

1 



A 



t3,(9,10) 



2 



g'^-7g^ + 18g-17±(g''-lV + 81g^-250g'^ + 442g^-436g + 193)^/^ . (2.21) 



Note that some Ajsj's vanish if g = 2 or g = 3; specifically, if g = 2, then Ata^s = and two 



of the roots of the cubic equation (|2 . 1 7| ) vanish, while for g = 3, At3_3 = and one of the 
roots of eq. (|2.17|) vanish. The vanishing of the chromatic polynomial at g = 0, 1, 2 and (i) 
g = 3 for the cyclic strip with m ^ mod 3 or (ii) g = 3 for the Mobius strip discussed 
below involves both some terms A^ vanishing and cancellations among others. This is also 
true of chromatic polynomials for strips of other lattices. 

Parenthetically, we note that the denominator of the generating function for this strip 
graph actually has dj) = 16: 

16 

V{trt{Ly = 3),FBCy,PBC^,q,x) = ^^(1 - A«j(g)x) (2.22) 

where the A^^- were given in eq. (|2.15|) and 

Ami = Am = -1 (2.23) 

Ai,i2 = Am3 = g - 3 (2.24) 
and Xtj, j = 14, 15, 16 are the roots of the cubic equation 

e + (2g - 5)(g - 3)^' + (g - 2fiq - 3)(g - 4)^ - (g - 2)^(g - 3) = . (2.25) 

A new feature of these results that we have not encountered before in our calculations of 
chromatic polynomials and generating functions for strip graphs is the occurrence of factors 
in T> that have multiplicity higher than unity: specifically, the factors (l + x)^, corresponding 
to \t3A = At3_ii = —1, and [1 — (g — 3)x]^, corresponding to Ai3^i2 = Xt3,i3 = g — 3. None of the 
At3,j, 12 < j < 16, contribute to P (and Xt3,u is identical to At3,i), so that P only involves 
the ten \t3,j's in eq. (|2.15|) . We had found cases earlier of some A/s in V not contributing 
to P, e.g., for the Ly = 2 cyclic strip of the kagome lattice and homeomorphic expansions of 



the Ly = 2 cyclic strip of the square lattice |11|. However, this is the first value of Ly where 



we have encountered this type of behavior for a strip of a homopolygonal lattice.[] 

^As in |3^, we define a homopolygonal 2D lattice as a regular tiling of the plane involving a single 
type of polygon, in which all vertices are equivalent, while a heteropolygonal lattice is the same but involving 
two or more different types of regular polygons. In the case of the kagome lattice, these are triangles and 
hexagons. 
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The numerator of the generating function has degree degx{M) = 14. The polynomials 
Aj, j = 1, ..,14 comprising this numerator can be calculated from the expressions for P and 
V, eqs. ( p.l5|) and (|2.22|) and hence are not listed here. 



From eq. ( |1.4|) , it follows that for m = 1 or 2 mod 3 and q = 3, the chromatic polynomial 
in eq. (|2.15|) vanishes; for m = mod 3 (and m > 3) and g = 3 we find 



P{tri{Ly = 3)m, FBCy, PBC^, g = 3) = 6 for m = mod 3 . (2.26) 

This may be contrasted with the case of a /c-critical graph Gk, where x{Cr) = k and P{G, q = 
x(G)) = k\ (e.g. /c = 3 for the triangular lattice). 
We define 

C{G) = CG, . (2.27) 

where the G-dependence in the coefficients is indicated explicitly. Note that for recursive 
graphs like the strip graphs considered here, the cgj depend on Ly and the boundary con- 
ditions, but not on L^. Our results above give 

C = q{q- 1)2 for G = tri{Ly = 3, FBCy, PBG^) . (2.28) 



in accord with the generalization |11, 12| C = P{Tiy,q) for this strip, where P{Tn,q) = 
q{q — 1)""^ and is the tree graph with n vertices. This was also true for the Ly = 3 cyclic 
strip of the square lattice lTl|, |l^ and is in accord with the coloring matrix approach [^, p7| . 

For the Ly = 3 triangular lattice Mobius strip (denoted by the subscript t3Mh) we 
calculate 

10 

P{tm{Ly = 3)„, FBGy, TPBG.,, q) = J2 Ct3MbA^t3Mb,jr (2-29) 

i=i 

where 

AfSMfej = Aj3j , j = I, .., 10 . (2.30) 

As was true in the case of the Lj^ = 2 Mobius strip of the triangular lattice, the coefficients 
Cj are not just of the c^'^^ polynomial form; instead, some coefficients are rational functions 
of q and others are algebraic nonpolynomial functions of q: 

Ct3Mb,i = c(') (2.31) 
= -^%(/_^^)^^ (2-32) 
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Ct3M6,(3,4) 



1 



(g-i)(g-5) 



Ct3Mb,5 



(g- -4g + 5) 
2(g-2) 



(2.33) 
(2.34) 



Ct3Mb,9 — QsMb.lO — C^^^ — 1 • (2.35) 

The j = 3,4 terms can be written in a manner which manifestly exhibits their property of 
being a symmetric function of Ajs 3 and Ai3^4: 



Ct3Mb,3{\3,3)"^ +Ct3Mb,4:{^t3,' 



(g-3)^ 
4fa-2) 



(At3,3)'" + (At3,4)'^ 



(g-l)(g-5)r(At3,3)™-(At3,4)' 



4(g - 2) 



At3,3 — Aj3,4 

(2.36) 

It is convenient to leave the expressions for Ct3Mb,j, j = 6,7,8, in the general forms that 
follow from eq. 



Ct3Mb,j — A 



t3,j 



^tm3,s(A 



t3,jj 



\9-s 



s=0 



n 



1 



l<i<10; i^j 



(At3j - At3,i) 



, J = 6,7,8 



(2.37) 



where the At3Mb,s are given in the appendix. 

For the sum of the coefficients for the Ly = 3 Mobius strip of the triangular lattice we 
calculate 

10 

C = J2ct3Mbj = q{q-l) for G = tri{Ly = 3, FBCy,TPBC,) (2.38) 

in agreement with a general formula that we have derived elsewhere for a Mobius strip graph 
Gs of width Ly of the square or triangular lattice |]35| , 



r for even Ly 

C{Gs{Ly), Mb) = I p{T^^^^ q) for odd Ly 



(2.39) 



Chromatic zeros for the Ly = 3, L^ = m = 2Q cyclic graph of the triangular lattice are 
shown in Fig. 0; with this value of m, the chromatic zeros for the triangular lattice Mobius 
graph are quite similar and both sets of chromatic zeros lie close to the boundary B and 
indicate its position. The locus E and the W functions are the same for the cyclic and 
Mobius graph families. We find 



q,{tTi{Ly = 3), FBCy, {T)PBC^) = 3 . 



(2.40) 
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Figure 2: B and chromatic zeros for the Ly = 3, Lx = m = 20 (i.e., n = 60) cyche triangular strip graph. 
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This will be discussed further below when we give our result for qc for the corresponding 
Ly = A strip of the triangular lattice. 

The locus B has support for Re{q) > and separates the q plane into five main regions. 
The outermost one, region extends to infinite |g| and includes the intervals q > 3 and 
g < on the real q axis. The other regions are labelled consequtively moving leftward along 
the real axis from q = qc = 3 to q = 0, and then including complex conjugate pairs of regions 
that have no overlap with the real axis. Thus, region R2 includes the real interval 2 < g < 3, 
while region includes the real interval < g < 2. There are also two complex-conjugate 
regions, R4, Rl centered at approximately q ^ 3.35 ± 0.7i. The triple points occur at (i) 
q ~ 3.14 + 0.37i and (ii) q ^ 3.12 + 1.12i, where -Ri, R2, and -R4 are contiguous, and (iii) 
q ^ 2.25 + 1.662, where R2, and i?3 are contiguous, together with the complex conjugates 
of these points. In addition, there are two tiny regions R^, R^ that touch the real axis at 
the single point qc = 3 and extend a very short distance to the upper and lower right, with 
triple points at about q ^ 3.02 ± 0.13i. These are plotted in Fig. ||. Thus, qc is a multiple 
point where four branches of B intersect. The occurrence of extremely small regions was also 
found for the {L^ 00 limit of the) Ly = 3 cyclic strip of the square lattice, as shown in 
Fig. 1 of |ll||. The fact that such tiny regions can occur means that when one performs the 
usual computer scan over the complex plane to map out the region diagram and determine 
the dominant Aj's in each region, since this involves a finite grid, one can only detect tiny 
sliver regions down to a certain resolution set by the grid of the scan. On the real axis, this is 
not a serious complication, since the mapping is a one-variable problem, but in the complex 
plane, one must a commensurately fine scanning grid to detect tiny sliver regions. 

The part of B and the associated chromatic zeros with largest magnitude occur at roughly 
q ~ 3.53 + 0.69i and have |g| ~ 3.59. Comparing the present locus B for Ly = 3 with the 
one obtained in [|l^ for Ly = 2, we obtain further evidence supporting the conjecture that 
for a strip graph of a regular lattice with periodic or twisted periodic longitudinal boundary 
conditions {B is the same for both of these), the envelope curve of B, i.e. the outermost 
portion of B for a given Ly surrounds the envelope curve for the B of the same strip with a 
smaller value of Ly. This was found to be true for the cyclic and Mobius strips of the square 
lattice for the cases Ly = 2,3 and 4 (and, for the cyclic case, also Ly = 1) 0, |lT], For 
strips of a regular lattice A, as Ly 00, one expects that the envelope curves defined for 
each Ly will approach a limiting curve, which is precisely the envelope for the boundary B 
of the full 2D lattice A defined via this limit. 

In region Ri, X^^ g is the dominant Xj (with appropriate choice of branch cut for evaluation 
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away from the negative real axis), so 

W = (A«,9)'/' , g G i?i . (2.41) 



The fact that this is the same as W for the (FBCj^,FBCa;) case fl^ is a general result 
In region i?2, ^^3,4 is dominant, so 

\W\ = \\t^A^'^ , g G i?2 (2.42) 

(in regions other than only \ W\ can be determined unambiguously [0). In region i?3, W 
is given by the largest (in magnitude) of the roots of the cubic ( p.l7|) , which we label Xa^'- 

\W\ = \\t^^^\^/\ qeRs. (2.43) 

In regions R4^,Rl, 

\W\ = \Xts,s\'^' , qeR,,Rl. (2.44) 

where A^s^g is the other among these cubic roots that has maximal magnitude in this region. 
Ini?5, Rl \W\ = \Xt3fi\'^\ 

3 Ly = A Cyclic Triangular Strip 

We have succeeded in performing an exact calculation of the chromatic polynomial for the 
cyclic triangular strip of arbitrarily great length and of the next larger width, namely, Ly = 4. 
The calculation is considerably more involved than for the Ly = 3 cyclic strip, as is indicated 
by the number of Aj terms in eq. (|1.3|) , namely, A^a = 26, as contrasted with the value Ly = 10 
for the Ly = 3 cyclic strip. Elsewhere we have proved that Nx is the same for the square 
and triangular strips of a given width and have given a general determination of A^a as a 
function of Ly [^. To go beyond this width for this family of cyclic strips of the triangular 



lattice would be increasingly difficult, since the number of terms A^a in (|1.3| ) is 70, 192, and 
534 for Ly = 5, 6, and 7, respectively. For = 4 we find the chromatic polynomial 

26 

P{tm{4 X m),FBCy,PBa,,q) = Q4,,(Am,,)"^ (3.1) 

where one term is Aj4,i = 1 and the others, A^j, are the roots of (i) the quartic equation 
( |9.2.1|) , for 2 < j < 5; (ii) the quartic equation ( |9.2.2|) , for 6 < j < 9; (iii) the degree-8 



equation ( |9.2.7|) , for 10 < j < 17; and (iv) the degree-9 equation (|9.2.16|) , for 18 < j < 26. 



Since these equations are somewhat lengthy, they are given in the appendix. In contrast to 
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the Ly = 2 and Ly = 3 cyclic triangular strips, here it is not possible to solve for all of the 
terms A^j as algebraic roots since for j = 10 through j = 26 these are roots of equations of 
degree higher than quartic. The corresponding coefficients cnj are 



Q4,i = c^'^ (3.2) 

ct^j = c(°) for 2 < j < 5 (3.3) 

c^^ . = c(3) for 6 < J < 9 (3.4) 

= c(2) for 10 < J < 17 (3.5) 

Cf4j = c^^) for 18 < J < 26 . (3.6) 



As before, one can, equivalently, list the generating function; however, since all of the relevant 
information is present in ( |3.1| ) with the requisite definitions of the coefficients Ct^j and terms 
At4j, we shall forego this. 

For the Mobius Ly = 4 strip of the triangular lattice, a general argument shows that 
the Aj's are the same as for the cyclic strip, although the coefficients cj are different and 
are more complicated than those for the cyclic strip, as was already shown by the Ly = 2 
triangular strip, where two coefficients for the Mobius case are algebraic nonpolynomial 
functions of q. As discussed before, since B is determined only by the Aj's, it is the same for 
the cyclic and Mobius strips of a given lattice. 

Chromatic zeros for the Ly = 4, L^ = m = IQ (hence n = 64) cyclic graph of the 
triangular lattice are shown in Fig. ^; the complex zeros give an approximate indication 
of B. The locus B separates the q plane into six main regions, including four that contain 
intervals of the real axis. The outermost one, region Ri, extends to infinite |g| and includes 
the intervals q > qc and g < on the real q axis. Region R2 includes the real interval 
3 < (? < region R3 includes the real interval 2 < g < 3, and region R4 includes the 
real interval < g < 2. There are also complex-conjugate regions -R5 and -R5 centered at 
approximately q ~ 3.5 ± 0.6i. There could also be other such complex-conjugate pairs of 
regions. The density of zeros is observed to be smaller on the parts of B extending through 
q = 2 and q = 3 and also on the right-most bubble-like curves than on the rest of the of B. 
As before with the Ly = 3 strip, it is straightforward to determine the various triple points. 
The maximal point at which the locus B crosses the real axis is 

qc{tri{Ly = 4), FBCy, {T)PBC^) - 3.228126 . (3.7) 

Thus, for the triangular strips with {FBCy, {T)PBCx) boundary conditions for which the 
chromatic polynomials have been calculated so far, i.e. for the widths Ly = 2, 3, and 4, qc 
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Figure 3: B and chromatic zeros for the Ly = 4, Lx = m = 16 (i.e., n = 64) cyche triangular strip graph. 
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is a nondecreasing function of Ly-. qc(tri{2)) = 3, qc(tri{3)) = 3, and qdtri^A)) = 3.228.... 
We have found that this is also true for the strips of the square lattice with the same 
cyclic or Mobius boundary conditions: q{sq{l)) = qc{sq{2)) = 2, qc{sq{3)) ^ 2.34, and 
qc{sq{4:)) ~ 2.49. This contrasts with the non-monotonic behavior of qc as a function of 
Ly that we have found for strips with {PBCy, FBCx) boundary conditions (see eq. ( ^.16|) 



below). It is anticipated that for cyclic and Mobius strips of the triangular lattice, as Ly 
increases beyond 3, qc will approach the 2D value qdtri) = 4. 

Another important feature of the locus B for the Ly = A cyclic triangular strip, which is 
the same as was true of the Ly = 2 and Ly = 3 cyclic triangular strips, is that it crosses the 
real axis at the points g = 0, 2 and 3. This leads us to the conjecture 

Conjecture : i3 D {g = 0, 2, 3} for tri{Ly, FBCy, {T)PBC^) V Lj, > 2 (3.8) 

The crossing at g = is present in a very wide class of families of graphs that contain global 



circuits p3[| . 



The crossing at g = 2 that we have found for Ly = 2,3,4 signals a zero-temperature 
critical point of the Ising antiferromagnet on these triangular strips, and similarly for the 



generalization ( |3.8|) to arbitrary Ly. Note that this T = critical point involves frustration, 
since not all of the spin-spin interactions around a triangle can have their energies minimized. 
Frustration is also responsible for the T = critical point of the Ising antiferromagnet on 



the full 2D triangular lattice but the nature of the critical singularities in the free 



energy, correlation length, etc. are different: they are algebraic for the triangular lattice, but 
are essential singularities for the Ly x oo strips . As will be shown below, this crossing of 
;B at g = 2 is also found for the Ly = 2 strips with torus or Klein bottle boundary conditions, 
i.e. {PBCy, PBCx) or {PBCy,TPBCx). Thus, it appears to be a general feature of the 
strip graphs of the triangular lattice with periodic or twisted periodic longitudinal boundary 
conditions. Indeed, the respective boundaries B for the Ly = 2 cyclic and Mobius strips of 
the square, triangular, and kagome lattices, and for the Ly = 3 cyclic and Mobius strips of 



the square lattice |]Tl| and the torus and Klein bottle strips of the square lattice pB] also 
cross the real g axis at g = 2 (as well as g = 0). 

In contrast, for strips with {FBCy, FBCx) (open) or {PBCy, FBCx) (cylindrical) bound- 
ary conditions, B does not pass through g = and, while B for the Ly = 3 open square strip 
crosses the real axis at g = 2, the respective loci B for the Ly = 4, 5 open and cylindrical 
square strips and the open and cylindrical triangular strips with Ly = 3,4,5 do not, as 
discussed further below. 

In regions Rj, 1 < j < 4, the dominant terms are, respectively, the root of maximal 
magnitude of (1) the first quartic equation (|9.2.1| ), the second quartic equation (|9.2.2|) , (3) 
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the eighth-degree equation ( |9.2.7|) , (4) the ninth-degree equation (|9.2.16|) . In region Ri, 



we label the dominant term as \t4,Ri, and so forth for the other regions. In the complex- 
conjugate regions R5, the dominant term is the largest-magnitude root of the ninth-degree 
equation (|9.2.16| ), which we label as Am,r5- We have 

W = (Ai4,m)'/' , q^Ri (3.9) 

\W\ = |Ai4,fl,f/^ , qeR, ,J = 2,3,4 (3.10) 

and 

\W\ = \Xu,n5\'^\ qeR,,Rl. (3.11) 

As usual, the boundaries of the regions are defined by the degeneracies, in magnitude, of the 
leading Xu/s in these regions. For example, on the real axis, one has |At4^/ji| = |Ai4,i?2| at 
Qc, |At4,_R2| = |At4,_R3| at g = 3, |At4_/j3| = |Ai4,R4| at g = 2, and |A(4,fi4| = |AM,ra| at g = 0. 

In contrast to the situation for the Ly = 2 and Ly = 3 cyclic strips of the triangular 
lattice, here Qc is a regular instead of a multiple (= intersection) point on the algebraic curve 
forming B, in the terminology of algebraic geometry (a multiple point on an algebraic curve 
is a point where several branches of the curve cross each other). 

Another qualitative difference is that for Ly = 2 and Ly = 3, B contained support only 
for Re{q) > 0, and the only point on B with Re{q) = was the origin itself; however, for 
Ly = 4, B extends slightly into the half-plane with Re{q) < 0. Although these lattice strips 
involve global circuits, it was found in that this is not a necessary condition for chromatic 
zeros and loci B of lattice strips (or homeomorphic expansions thereof) to have support for 



Re{q) < 0; this was also evident in results of ||29[, ruling out a conjecture in ^7j. In the case 
of strips of the square lattice we had found that for Ly = 1 and 2, B only had support for 
Re{q) > 0, but for Ly = 3 and Ly = 4 it had support also for Re{q) < 0. 

Finally, another general feature is that the motion of the outermost curves forming part 
of B for Ly = 2, 3, 4 is consistent with the inference that as Ly increases, these form a limiting 
curve. 

4 Ly = 3 Triangular Lattice Strips with Torus and Klein 
Bottle Boundary Conditions 

For the Ly = 3 triangular lattice strips with torus boundary conditions (denoted with the 
subscript tt), we calculate, by the same methods as above, 

11 

P{tri{Ly = 3U PBCy, PBC, q) = Y, c«j(A«,,)"^ (4.1) 
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where 



Xu,2 = q-2, Cit,2 = ^(g- i)(g-2) , 



Att,3 = 3g - 14 , cu,3 = ^q{q - 3) , 



Xu,^ = -2{q - Af 



cttA = g - 1 



(4.2) 
(4.3) 
(4.4) 
(4.5) 



Ah,5 = - 9g' + 29g - 32 



Qi,5 = 1 



(4.6) 



Xtt,{6,7) 



^±27ri/3 



Ctt,6 — Qt,7 



-g(g-l)(2g-7) 

D 



(4.7) 



Att,(8,9) = (g^ - 5g + 7)e 



±27rj/3 



(4. 



Att,(io,ii)) = -(2g - 7)e 



±27rj/3 



Cu,io = Cu,n = -iq-l)iq-2) . (4.9) 



Note that ahhough several of the coefficients are complex (and have coefficients that are not 
integers), the chromatic polynomial itself is, of course, a polynomial with integer coefficients. 
The sum of the complex terms can be written as 



11 



j=6 



q{2q - 7) 



f2(g2-5g + 7)'" + (g-2)[-(2g-7)]' 



f 217171 



(4.10) 



None of these terms with 6 < j < 11 can be dominant anywhere in the q plane since this 



would imply that the limit ( |1 . 21 ) would not exist. As compared with previously calculated 
strip graphs, this exhibits a number of qualitatively new features: (i) previously, only one Xj 
was a constant, independent of q, and this was always equal to either 1 or —1; here there are 
three, and none of these is equal to ±1; (ii) previously, all of the A^'s were either polynomials 
with real coefficients or algebraic functions of polynomials with real coefficients. By general 
arguments given before ||23[, the dominant Xj in region Ri has as its highest power g"/*"^ i.e., 
in the present case, g^, and has coefficient 1. We find that the generating function for this 
case has degxi^hf) = 9 and degxiV) = 11, so that all of the A/s in V contribute to P. 
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Given that x = 4 for m = 1 or 2 mod 3 from eq. (|1.5|) , it follows that for these m values, 
and q = 3, P{tri{Ly = 3)^, PBCy, PBC^, g = 3) = 0; however, x = 3 for m = mod 3; in 
this case, we find (for m > 3) 

P{tri{Ly = 3)^, PBCy, PBC^, g = 3) = 6 for m = mod 3 . (4.11) 



This is analogous to eq. (|2.26| ) above 



For the Ly = 3 triangular lattice strips with Klein bottle topology (denoted with the 
subscript tk) we calculate 



5 

P{tri{Ly = 3)m,PBCy,TPBC^,q) = ^ Qfcj(Atfc,,)"^ 



= -(g - l)(-2r - ^(g - l)(g - 2)(g - 2)™ + ^g(g - 3)(3g - 14)™+ 

(g - 1) [-2(g - 4)2]™ + (g^ - 9g2 + 29g - 32)™ . (4.12) 

Thus, Ctk,j = Ctt,j, j = 1, ..S, and the six Xuj, j = 6, ..11 that involve complex factors do not 
contribute to P(tri{Ly = 3)m, PBCy, TPBCx, g). Since none of these six \u/s is dominant 
anywhere, it follows that B is the same for the Ly = 3 triangular lattice strips with torus 



and Klein bottle topology, just as was true of the analogous square lattice strips ||2^, and 



in accord with the general discussion of Ref. |2^. For this case the generating function has 
degx{M) = 2 and degx{T>) = 5, so that, as was true of the Ly = 3 Mobius and torus strips 
of the triangular lattice, all of the A/s in V contribute to P. 

For the sum of the coefficients we find, for the triangular Ly = 3 strip with torus boundary 
conditions 

C = P{K3, q) = q{q - l)(g - 2) for G = trt{Ly = 3) with {PBCy, PBC^) (4.13) 

and for this strip with Klein bottle boundary conditions 

^ = for G = tri{Ly = 3) with {PBCy,TPBCx) ■ (4.14) 

These are the same respective values as those calculated in for the corresponding Ly = 3 
strips of the square lattice with torus and Klein bottle boundary conditions. 

Chromatic zeros for the Ly = 3, L^ = m = 20 torus graph of the triangular lattice are 
shown in Fig. ^ as before, with this value of m, the chromatic zeros for the torus and Klein 
bottle boundary conditions are generally similar and both sets of chromatic zeros lie close 



20 



4.0 



3.0 



2.0 



1.0 



im(q) 0.0 



-1.0 



-2.0 



-3.0 



-4.0 I 
-1.0 



0.0 



1.0 



2.0 
Re(q) 



3.0 



4.0 



5.0 



Figure 4: B and chromatic zeros for the Ly — 3, Lx = m = 20 (i.e., n = 60) triangular strip graph with 
torus boundary conditions. 
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to the boundary B and indicate its position. The locus B and the W functions are the same 
for the torus and Klein bottle graph families. We find 



Qc ~ 3.7240756... for {G} = tri{Ly = 3) with {PBCy, PBC^) or {PBCy, TPBC^) . 

(4.15) 

This value of Qc is given by the real root of the degeneracy equation — 9q^ + 32g — 46 = 0. 
It is interesting that this is just 7 % below the value for the infinite 2D triangular lattice, 
qdtri) = 4. Thus, for the width Ly = 3 strips where a comparison can be made, the use of 
both periodic transverse boundary conditions and periodic or twisted periodic longitudinal 
boundary conditions expedites the approach to the 2D thermodynamic limit, in the sense 
that it yields a value of qc that is substantially closer to the 2D value than was obtained 
for the same width strip with free transverse boundary conditions in eq. ( |2.40| ). This is 



expected since for the torus or Klein bottle boundary conditions (a) the resultant graphs are 
A-regulaiQ with the degree (coordination number) A = 6 of the 2D triangular lattice, and 
(b) there are no boundaries to the surface on which the graphs are embedded. In contrast, 
for the boundary conditions of types (i)-(iv) (aside from degenerate cases) the graphs are 
not A-regular and they do have boundaries. 

The locus B has support for Re{q) > and separates the q plane into three regions. The 
outermost one, region Ri, extends to infinite \q\ and includes the intervals q > qc and g < 
on the real q axis. Region R2 includes the real interval 2 < q < qc and extends upward and 
downward to the complex conjugate triple points on B at qt,qt — 4.0 ± 1.7i. Region R^ is 
the innermost one and includes the real interval < g < 2. The maximum value of \Im{q)\ 
on the boundary between Ri and -R3 is about 3.0, occurring at Re{q) ~ 2.6. The maximum 
value of |g| on B is roughly 4.6, occurring at g ~ 3.9 ± 2.5i. The boundary between R2 
and i?3 curves to the right as one increases \Im{q)\, extending from q = 2 upward to qt and 
downward to q^. As is evident in Fig. ^, the density of chromatic zeros is highest on the 
Ri — R2 boundary near qc- The feature that B has support only for Re{q) > is the same 
as was found p6[ for the Ly = 3 strips of the square lattice with torus and Klein bottle 



topology. 

In region Ri, \tt,5 = Xtk,5 is the dominant Xj, so 

W = {q^ - 9g2 + 29g - 32)^/^ , q E Ri . (4.16) 



The fact that this is the same as W for the (PBCy,FBC2;) case, eq. (|5.3| ), is a general result 
23| . The importance of the PBCy is evident from the fact that for the same width of three 



squares, the strip with (FBCj,,FBCa;) yields a different W [|T7 



A-regular graph is a graph all of whose vertices have the same degree, A, where the degree of a vertex 
is defined as the number of edges connected to it. 
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In region R2 Xtt,3 = Xtk,3 is dominant, so 

\W\ = \3q - U\^/^ , qeR2 (4.17) 

(in regions other than only \W\ can be determined unambiguously 0). In region i?3, 
Att,4 = Atfc,4 is dominant, so 

= |2(g- 4)^/3 , gei?3. (4.18) 
At all of the three points, g = 0, 2, and q = qc — 3.72 where B crosses the real q axis, it does 



so vertically. The present results are in accord with the inference [|T7|, [TT[ that for a recursive 
graph with regular lattice structure, a sufficient condition for B to separate the q plane into 
two or more regions is that it contains a global circuit, i.e. a path along a lattice direction 
whose length goes to infinity as n — > 00; here this is equivalent to PBC^. 

Our calculations of the zero-temperature Potts antiferromagnet partition functions (chro- 
matic polynomials) and exponential of the entropy, W, for Ly = 3 strips of the triangular 
lattice with various boundary conditions including free and periodic or reversed-orientation 
periodic thus elucidate the role that these boundary conditions and the associated topologies 
play. One particular feature of note is that the torus and Klein bottle graphs have interest- 
ingly different chromatic polynomials, with different A^a, but the W functions and hence the 
boundaries B are the same. 

5 Width = 5, 6 Strips of the Triangular Lattice with 

{PBCy, FBC,^ 



We first briefiy review the cases Ly = 3,4 |T8|, ^ |2^. The chromatic polynomial for the 
Ly = 3 triangular (t) strip with n = 3(m + 2) vertices (following the labelling conventions in 
p!7| ) and {PBCy, FBC^) boundary conditions, denoted t3PF in the subscripts, has Nx = 1, 

At3PF = - + 29g - 32 , (5.1) 

P{trt{3 X m, PBCy, FBC^), q) = q{q -l){q- 2){Xt3PFr^' , (5.2) 



and 



with i3 = 0. 



W{tri{Ly = 3), PBCy, FBC^, q) = {q^ - 9q^ + 29q - 32)^/^ (5.3) 
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For the Ly = 4 strip with [PBCy, FBCx) boundary conditions and n = 4(m + 2) vertices, 



one had [||, ^ Nx = 2 and 



where 



where 



P{tri{A X m, PBCy, FBC,), q) = J2 CmpfjIAmpfJ'"^' (5.4) 



X _(?-3) 

MiPFJ — 



2 



(5.5) 



Tj4p^ = - 9g2 + 33g - 48 (5.6) 
R^^pp = q^- lOq^ + 43g^ - 106g + 129 . (5.7) 
In this case, as is evident in Fig. 4(a) of |13]), B includes both arcs and a closed oval invariant 



under complex conjugation which crossed the real axis where the prefactor of the square root 
vanished, i.e. at Qc = 4, and at the real root of Tupfi namely q ~ 3.481406. There are two 
regions, Ri and R21 the exterior and interior of the oval. We have 

W{tri,A X 00, PBCy, q) = {Xtmf^ , 9 ^ i?i (5.8) 

and 

W{tri, 4 X 00, PBCy, q) = (Ampp,2)'/^ , q e R2 (5.9) 

with appropriate choices of branch cuts. 

We have calculated new results for the width Ly = 5 and Ly = 6 strips with {PBCy, FBC^) 
boundary conditions and n = Ly{m + 2) vertices. For Ly = 5 we find A^a = 2 and 

2 

P{trt{5 X m, PBCy, FBC,), g) = ^ Q5PPj(Ai5PP,,)™+' (5.10) 



where 



where 



At5PF,j - 2 



and 



Tt5PF ± St5PF \/ RtbPF (5.11) 

Tt5PF = - 15g^ + 98g^ - 355g^ + 711g - 614 (5.12) 

R^^pp = q^- 8q^ + 26g^ - 60g + 85 (5.13) 

^t5PP = g^-llg' + 43g-58 . (5.14) 
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As before, the coefficients Ct^pFj can be computed using eq. (|2.9| ) in terms of the generating 
function. This generating function is given in the appendix. 

Just as was true for the Ly = 4 triangular lattice strip with {PBCy, FBCx), in the present 
Ly = 5 case B includes both a pair of complex conjugate arcs and an oval that crosses the 
real axis at two points. Analogously to the previous case, one of these points is the real zero 
of Tt5PF at 

Qi ~ 3.207224 (5.15) 
and the other is the real zero of the pref actor St5PF in front of the square root, at 

q = 3.284775 = for {G} = tri{5 x oo, PBCy, FBC^) . (5.16) 

Taking the approximate center of the oval as the average of the left and right crossings, 

Qcenter = (?£ + Qc) , OUe fiuds that 

g,enter ^ 3.2460 (5.17) 
near to the Tutte-Beraha number B^ ^ 3.246980.., where [^, [14 



Br = 2 + 2cosf— ) . (5.18) 



\ r J 



Chromatic zeros near Bj were also noted in [g^ for these types of strips. Note that q^ 
decreases from 4 to the above value in eq. ( ^.16 ) as Ly increases from 4 to 5. Hence, in con- 
trast to our calculations of strips of the square and triangular lattices with {FBCy, FBCx), 
{FBCy, (T)PBCx), and {PBCy, (T)PBCx), where in all cases considered, the respective 
gc's were nondecreasing functions of Ly, here we find that this is not the case for boundary 
conditions of the type {PBCy, FBCx). 

Chromatic zeros and B for this strip are shown in Fig. ^; B separates the q plane into two 
regions, Ri, the exterior, and R2, the interior, of the oval shown in the figure. Comparing 
the ovals for the Ly = 4 strip (Fig. 4(a) of []T^) and the Ly = 5 strip (Fig. |^ here), one sees 
that the oval shrinks in size. We have 

W{tri, 5 X 00, PBCy, q) = (At5PF,i)'/' , q e Ri (5.19) 

and 

W{tm, 5 X 00, PBCy, q) = (At5PF,2)'/' , q e R2 (5.20) 
with appropriate choices of branch cuts. 

For the strip of with Ly = 6 and {PBCy, FBCx) we find Nx = 5 and 

5 

P{trz{Q X m,PBCy,FBCx),q) = Q6PPj(Ai6PF,,)"^+' (5.21) 
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Figure 5: Locus B for the width Ly = 5 strip (tube) of the triangular lattice with {PBCy, FBCx) boundary 
conditions. Thus, the cross sections of the tube form pentagons. For comparison, chromatic zeros calculated 
for the strip length = m + 2 = 16 (i.e., n = 80 vertices) are shown. 
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Here the A(6pf's are the roots of a fifth-order equation, so we cannot solve for them in terms 
of algebraic expressions as was possible for previous strips of this type with Ly = 3, 4, 5. We 
give the generating function for this strip in the appendix. 




1 2 3 4 5 



Re(q) 



Figure 6: Locus B for the width Ly = 6 strip (tube) of the triangular lattice with {PBCy, FBCx) boundary 
conditions. Thus, the cross sections of the tube form hexagons. For comparison, chromatic zeros calculated 
for the strip length ^ m + 2 — 16 (i.e., n = 96 vertices) are shown. 

Chromatic zeros for this strip are shown in Fig. |^ for = m + 2 = 16, i.e., n = 96; this 
length is sufficiently great that the chromatic zeros give an approximate indication of the 
location of the locus B. The locus thereby inferred appears to consist of a single connected 
set of curves and crosses the real axis. From our exact analytic expressions, we calculate the 
this crossing point to be 

= 3.252419 for {G} = tri{6 x oo, PBCy, FBC^) . (5.22) 

The morphology of chromatic zeros for this long 6 x 16 cylindrical strip is similar to that 
found for a 8 x 8 patch with cylindrical boundary conditions in p9[. In both cases, one 



can figuratively think of how the inferred loci B can be modified to yield the locus for the 
Ly = oo strip with {PBCy, FBCx) cylindrical boundary conditions: one (i) pulls the left- 
hand endpoints further to the left (into the Re{q) < half-plane) and around so that they 
meet in a cusp at the origin, q = 0; (ii) pulls the right-hand endpoints of the prongs over 
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and around so that they meet in a cusp at g = 4; (iii) shifts the crossing curve shghtly to the 
right so that it crosses the real axis at g ^ 3.82; and (iv) makes minor shifts of the rest of the 
locus B so as to obtain B for the cylindrical strip with Ly = oo. In this respect, just as was 



discussed in , one sees that the periodic transverse boundary conditions help to minimize 
finite-size effects. Clearly, for a fixed Ly, with L^^ ^ oo, finite-size effects are reduced the 
most for (PBCy, {T)PBC^), next-most for {PBCy, FBC^) and {FBCy, {T)PBC^), and least 
for {FBCy,FBC^). We find 

W {tri, 6 X oo, PBCy, q) = {XtePF,maxY^'' for q e Ri (5.23) 

where Xt6PF,m.ax denotes the root with the maximal magnitude. 

It is of interest to use these results to study the approach of W to the limit for the full 
infinite 2D triangular lattice, extending the work of |2y]. In Table |1] we list the various values 
of W{tri{Ly X oo, PBCy, FBCx),q), denoted as W{Ly,q) to save space, together with the 
corresponding values of W for the full 2D triangular lattice, W{tri, q), obtained via numerical 



evaluation of an infinite product representation in [2£] as checked by series expansions and 



rigorous bounds [^, and the ratio Rw{tri{Ly x oo), PBCy, FBCx,q), where 



^ / ./v X ^ W(tri(Ly X oo),BCy,BCx,q) 

Rw{tri{Ly X oo), BCy, BC,, q) = • ^ • (5-24) 

This extends the previous study in Recall that for g > 4, W{tri{Ly x oo), BCy, BC^, q) 
is independent of BCx pO| , p3| . Evidently, the approach of W for the infinite- length finite- 
width triangular strips to the 2D thermodynamic limit as Ly increases is quite rapid; for 
moderate values of q, say 6 or 7, the ratio Rw(tri{5 x oo), PBCy, FBC^, q) is equal to 1 to 
within approximately 10~^ or better, and the ratio Rw(tri{6 x oo), PBCy, FBC^.q) is at 



least as close to 1. As was proved in the approach is non-monotonic for PBCy (and 



monotonic for FBCy) 



6 Width Ly = 5 Strips of the Triangular Lattices with 

{FBCy, FBC) 

Previous studies have been published on strips of various lattices with {FBCy, FBCx) bound- 
ary conditions [1^, |19|. In the case of the the square and triangular lattice, these went 
up to = 4 and involved A/s that were roots of a cubic and quartic equation, respec- 
tively. Although one can calculate chromatic polynomials for wider strips, the analysis is 
more cumbersome if the equations defining the A/s are higher than quartic, so that one 
cannot solve for these Aj's as analytic closed-form algebraic expressions. A study of these 
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Table 1: Comparison of values of W{tri{Ly), PBCy,q) with W{tri, q) for 4 < g < 10 and BCx = FBCx or 
{T)PBCx- For each value of q, the quantities in the upper line are identified at the top and the quantities 
in the lower line are the values of Rw{tri{Ly), PBCy, q). The PBCy is symbolized as Py in the table. 



Q 


W{trt{3),Py,q) 


W{tm{4),Py,q) 


W{tm{5),Py,q) 


W{tm{6),Py,q) 


W{tri, q) 


4 


1.58740 
1.0865 


1.18921 
0.8140 


1.39252 
0.9531 


1.49603 
1.0240 


1.46100 
1 


5 


2.35133 
1.0385 


2.21336 
0.9776 


2.26877 
1.0021 


2.26894 
1.0021 


2.26411 
1 


6 


3.23961 
1.0112 


3.185055 
0.9941 


3.20718 
1.0010 


3.20399 
1.0000 


3.20388 
1 


7 


4.17934 
1.0027 


4.15965 
0.99795 


4.16987 
1.0004 


4.16805 
1.0000 


4.16819 
1 


8 


5.14256 
0.99980 


5.13936 
0.9992 


5.14446 
1.0002 


5.14348 
1.0000 


5.14358 
1 


9 


6.11803 
0.99880 


6.12324 
0.99965 


6.12587 
1.0001 


6.12533 
1.0000 


6.12539 
1 


10 


7.10059 
0.99849 


7.11027 
0.99985 


7.11161 
1.0000 


7.11131 
1.0000 


7.11134 
1 



on wider open strips is in . As an illustration of this sort of situation, here we present a 
calculation of chromatic polynomials for the strip of the triangular lattice with width Ly = 5 
and {FBCy, FBCx). For this width, the A/s are solutions of a degree-9 equation and hence 
cannot be solved for as algebraic roots. The method of calculation is again the iterated 
use of the deletion-contraction theorem. In |l^, a given strip {Gs)m was constructed by m 



successive additions of a subgraph H to an endgraph J; here, I = H, so that, following the 



notation of |T^, the total length of the strip graph {Gs)m is m + 2 vertices, or equivalently, 
m + 1 edges in the longitudinal direction. The results are conveniently expressed in terms 
of the coefficient functions in the generating function, as discussed above. For the width 
Ly = 5 open strip of the triangular lattice we find degzi'D) = Nx = 9. The coefficient 
functions h^FFj and Ai^ffj (cf. eqs. ( ^^) and (|2.3| )) are listed in the Appendix, where FF 



is short for {FBCy,FBCx). 

In Fig. we show a plot of chromatic zeros for the open strip of the triangular lattice 
with Ly = 5 and length m + 2 = 16 vertices, so that the strip has n = 80 vertices in all. 
From an analysis of the degeneracy of leading Aj's, we find the exact result 

qc = 3 for sq{5x oo,FBCy,FBCx) . (6.1) 
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Figure 7: Chromatic zeros for the Ly = 5 open strip of the triangular lattice of length m + 2 = 16 vertices 
(i.e. total number of vertices n = 80). 
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This is in agreement with the chromatic zeros for the finite strip shown in Fig. |^. Again, 
we can compare Fig. |^ with the corresponding plots for Ly = 2 and Ly = 3 (Fig. 5(a,b) of 



T7|), and, just as was true for the corresponding three open strips of the square lattice, this 
comparison shows that as Ly increases, the arcs forming B elongate and the arc endpoints 
nearest to the origin approach more closely to the origin. We recall that for the Ly = 4 open 
triangular strip, no qc is defined since B does not cross the real axis. 

For q > qc, we have, for the physical ground state degeneracy per site of the g-state Potts 
ant if erromagnet , 

W{t{5 X OO, FBCx, FBCy), q) = {\tr^,,,ma.f''^ (6.2) 

where Xgqj^max denotes the solution of eq. ( |2.8| ) with the coefficients ( |9.4.1| )-( p.4.9| ) that has 
the maximal magnitude in region Ri. 

As before, it is of interest to use this result to study the approach of W to the limit for 
the full infinite 2D triangular lattice, extending the work of for this set of boundary 
conditions. In Table ^ we list the various values of W{tri{Ly x oo, FBCy, BCx), q), denoted 
as W{Ly, q) to save space, together with the corresponding values of W for the full 2D 
triangular lattice, W{tri,q) and the ratio Rw(tri{Ly x oo), FBCy, BCx, q) = W(tri{Ly x 
oo), FBCy, BCx, q)/W(tri, q). Recall that the value of for g > 4 is independent of BCx 
pO| . Again, the approach of W for the infinite-length finite- width triangular strips to the 
2D thermodynamic limit as Ly increases is quite rapid. 



7 Comparative Discussion 

In this section we give a general discussion of the locus B. We have found several interesting 
features: 

1. For the strips of the triangular lattice of width Ly = 3 and Ly = A with {FBCy, {T)PBCx) 
and strips of width Ly = 3 and {PBCy, (T)PBCx) boundary conditions studied here, 
we have shown that the locus B encloses regions of the q plane including certain inter- 
vals on the real axis and passes through g = and 2 as well as other possible points, 
depending on the family. This extends the previous study of the Ly = 2 strip of the tri- 
angular lattice with {FBCy, {T)PBCx) While the strips with {FBCy, {T)PBCx) 
have a locus B that passes through q = 3, this is not the case with at least the Ly = 3 
strip with {PBCy, {T)PBCx). For the strips with {FBCy, {T)PBCx) considered here 
we have shown that q^ is a nondecreasing function of Ly. As one increases Ly, the 
Ly = 4 cyclic strip is the first one for which one can no longer solve for all of the Aj's 
as algebraic roots, since some of the equations involved are higher than quartic. 
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Table 2: Comparison of values of W{tri{Ly x oo),g) for FBCy and BCx = FBCx or {T)PBCx with 
W{tri,q) = W{tri{oc- x oo),(jf) for 3 < g < 10. For brevity of notation, we omit the {FBCy, BC^) in the 
notation. For each value of q, the quantities in the upper line are identified at the top and the quantities in 
the lower line are the values of Rw{tri{Ly), q). 



Q 


W{tri{2),q) 


W{tri{3),q) 


W{tri{'i),q) 


W{tri{5),q) 


W{tri{oo), q) 




4 


2 


1.77173 


1.67619 


1.62270 


1.46100 






1.369 


1.213 


1.147 


1.111 


1 




5 


3 


2.72998 


2.60495 


2.53251 


2.26411 






1.325 


1.206 


1.151 


1.1185 


1 




6 


4 


3.71457 


3.579715 


3.50112 


3.20388 






1.248 


1.159 


1.117 


1.093 


1 




7 


5 


4.70571 


4.56515 


4.48283 


4.16819 






1.200 


1.129 


1.095 


1.075 


1 




8 


6 


5.69974 


5.55530 


5.47040 


5.14358 






1.167 


1.108 


1.080 


1.0635 


1 




9 


7 


6.695395 


6.54810 


6.46129 


6.12539 






1.113 


1.093 


1.069 


1.055 


1 




10 


8 


7.69208 


7.54259 


7.45430 


7.11134 






1.125 


1.082 


1.061 


1.048 


1 
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2. The crossing of B at the point q = 2 for the (infinite-length hmit of) strips with global 
circuits nicely signals the property that the Ising antiferromagnet has a frustrated 



zero-temperature critical point on these strips. This has been discussed in in the 
context of exact solutions for finite-temperature Potts model partition functions on 
the Ly = 2 cyclic and Mobius triangular strips. In contrast, this connection is not, in 
general, present for strips with free longitudinal boundary conditions since B does not 
pass through g = 2 (as recapitulated below). The g = 3 Potts antiferromagnet also has 
a zero-temperature critical point on these strips, and this is similarly manifested by 
the crossing of the singular locus B through the point g = 3 for the cyclic and Mobius 
strips with {FBCy, (T)BCx), but not for the Ly = 3 torus and Klein bottle strips with 
iPBCy,iT)PBC,). 

An interesting feature of the cyclic strips of the triangular lattice (as well as the square 
and kagome lattices) is that for all of cases that have been studied here and in 0, |ll], |12 



there is a correlation between the coefficient cq^j of the respective dominant Ag^j's in 
regions that include intervals of the real axis. Before, it was shown ^3) that the cg,j 
of the dominant XgsJ region Ri including the real intervals q > qc{{G}) and g < is 
c^^^ = 1, where the c*^°'^ were given in eqs. (|2.10| ), ( |2.11|) . Here we extend this, observing 



that the cg^j that multiplies the dominant XgsJ the region containing the intervals 
< g < 2 is c^^\ For the cyclic strips of the triangular lattice that we have studied, 
namely, Ly = 2,3,4, the cg^j multiplying the dominant Xg^j in the region containing 
the interval 2 < g < 3 is c*-^-*. For the Ly = A strip of the triangular lattice, there is 
another region containing the real interval 3 < g < gc, where gc for this strip was given 
in (|3.71 ), and we find that c'-^^ multiplies the dominant A^j containing this interval. 
For the cyclic strips of the square lattice, although the values of gc are different, a 
similar correlation is observed; in particular, for the respective widths Ly = 3,4, c^^-* 
multiplies the XgsJ that is dominant in the region containing the interval 2 < q < qc, 



where gc — 2.34 and gc — 2.49 for these two widths 



4. For the strip of the triangular lattice of width Ly = 5 with cylindrical {PBCy, FBC^) 
boundary conditions, we find that B consists of arcs together with a (closed) oval. 
However, B does not pass through g = 0, 2, or 3. This is qualitatively the same 
morphology that was found for the corresponding strip of width Ly = 4 (see Fig. 4 
of ||18|). For comparison, the triangular strip with width Ly = 3 and {PBCy, FBCx) 
had i3 = 0. The point gc = 4 for Ly = 4 and gc — 3.28 for Ly = 5, which shows 
that for this family of strips, in cases where there is a gc (there is none for Ly = 3), 
it is not, in general, a nondecreasing function of Ly. This is somewhat reminiscent of 
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the non-monotonicity of W that we showed in the case of [PBCy, BCx) strips, where 
BC^ = FBC^ or (T)PBC^, in |g. For the Ly = Q strip with {PBCy, FBC^), we infer 



from the chromatic zeros that B consists of a single component, and from the analytic 
results we compute that qc — 3.25, which again shows the nonmonotonicity of Qc as a 
function of Ly. The morphology of chromatic zeros for our 6 x 16 strip is similar to 



that found in p9[ for 8x8 patch, both with cylindrical boundary conditions. For strips 
with torus or Klein bottle boundary conditions, W and B have been calculated for only 
one Ly value, namely, Ly = 3 for the square lattice in and for the triangular lattice 
here. For the Ly = 3 {PBCy, {T)PBCx) square strip, qc is equal to the value 3 for 
the full 2D square lattice, but for the Ly = 3 {PBCy, {T)PBCx) triangular strip, our 



result (|4.15|) shows that qc is less than the value of 4 for the 2D triangular lattice. 



5. We have included an illustrative result for a wider strip with {FBCy, FBCx) boundary 
conditions, namely the strip with Ly = 5. As one increases Ly, this is the first value 
at which one can no longer solve analytically for the Aj's as algebraic roots. The locus 
B is similar to the respective loci that were found earlier in studies of open strips 



1^, |T8|, |T9[ in that it does not pass through g = or g = 2 and does not separate 
the q plane into regions containing intervals of the real axis. (It does pass through 
g = 3, unlike the loci for the open strips with Ly = 3,4.) Our results confirm the 
trend that was observed in earlier work [l^, namely that as Ly increases, the arcs 
elongate and move closer together, and the arc endpoints nearest to the origin move 
toward this point. In contrast to the strips with global circuits, these strips do not 
manifest the property that the q = 2 (Ising) and q = 3 Potts antiferromagnets have 
zero-temperature critical points since B does not pass through these respective points. 
The simplest example of this is the Ising antiferromagnet on the infinite open line; in 
this case, although the model has a well-known zero-temperature critical point, this is 
not evident in the locus B, which is the emptyset. 

6. For the strips with {PBCy, FBC^) and {FBCy, FBCx), while the endpoints of the 
arcs on B that lie closest to the origin tend to move toward the origin as Ly increases, 
leading one to expect that in the limit — > oo, the limiting locus B would pass 
through g = 0, no such motion toward g = 2 is observed in the cases so far calculated. 
This is in agreement with the fact that the locus found in for the triangular lattice 
constructed as the limit L^, Ly ^ oo with {PBCy, FBCx) boundary conditions passes 
through g = and 4 (and at g ~ 3.82) but not through g = 2 or g = 3. Thus, assuming 
that our conjecture ( p.8|) is correct, it follows that the locus B for the triangular lattice 
depends on the boundary conditions used to define this lattice: if one constructs as 
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the limit Lj, — > oo with {PBCy, FBCj) (cyhndrical) boundary conditions, then B does 
not pass through g = 2 or 3 p9|, while if one constructs it as the limit Ly ^ oo with 
{FBCy, (T)PBCx) (cyclic or Mobius) boundary conditions, then, if the conjecture is 
valid, B would pass through q = 2 and 3 in the limit just as it does for each of the 
values Ly = 2,3,4 studied so far. However, since the value of qc pertains directly to 
a physical quantity, as the minimal value of real q above which W{q) is analytic, one 
expects that qc should be independent of the boundary conditions used to define the 
2D lattice. All results obtained so far are consistent with this expectation. 

7. There have been a number of theorems proved concerning real chromatic zeros. An 
elementary result is that no chromatic zeros can lie on the negative real axis g < 0, 
since a chromatic polynomial has alternating coefficients. It has also been proved that 
there are no chromatic zeros in the intervals < g < 1, and 1 < g < 32/27 [Q. The 
bound of 32/27 in |^ has been shown to be sharp; i.e., for any e > 0, there exists 



a graph with a chromatic zero at g = 32/27 + e [S]. See also [^, ^ and references 
therein. Based on our studies of strips of the triangular (and square) lattices with the 
various boundary conditions considered, we make the following observation: for such 
strips, we have not found any chromatic zeros, except for the zero at g = 1, in the 
interior of the disk |g — 1| = 1. This motivates the conjecture that for these strips, 
there are no chromatic zeros with |g — 1| < 1 except for the zero at g = 1. Assuming 
that this conjecture is valid, the bound would be a sharp bound, since the circuit graph 
with n vertices, C„, has chromatic zeros lying precisely on the circle |g — 1| = 1 and at 
g = 1 [|^. Further work is needed to prove (or disprove) this conjecture. 

Some relevant features are summarized in Table |^. The entries for Ly = oo with 



{PBCy, FBCa) are from P . 



8 Conclusions 

In this paper we have presented exact calculations of the zero-temperature g-state Potts 
antiferromagnet partition functions (equivalently, chromatic polynomials P), on strips of 
the triangular lattice of width Ly = 3 and with boundary conditions of four types: (a) 
{FBCy, PBC^) = cyclic, (b) {FBCy, TPBC^) = M5bius, (c) {PBCy, PBC^) = toroidal, and 
(d) {PBCy, TPBCx) = Klein bottle, where F, P, and TP denote free, periodic, and twisted 
periodic. In the infinite-length limits of these strips, exact results were given for the ground 
state degeneracy (exponential of the ground state entropy), W , and its analytic structure in 
the complex g plane, in particular, the nonanalytic locus was discussed. Exact calculations 
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Table 3: Comparative listing of properties of chromatic polynomials P, groimd state degeneracy functions 
W, and nonanalytic loci B for strip graphs Gg of the triangular (tri) lattice and their infinite-length limits. 
New resiilts in this work are marked with asterisks in the first column. The properties apply for a given 
strip of type Gg of size Ly x L^; some apply for arbitrary L^, such as N\, while others apply for the 
infinite-length limit, such as the properties of the locus B. For the boundary conditions in the y and x 
directions {BCy, BC^), F, P, and T denote free, periodic, and orientation-reversed (twisted) periodic, and 
the notation (T)P means that the results apply for either periodic or orientation-reversed periodic. The 
column denoted eqs. describes the numbers and degrees of the algebraic equations giving the XgsJ in P't for 
example, {3(1), 2(2), 1(3)} indicates that there are 3 linear equations, 2 quadratic equations and one cubic 
equation. The column denoted BCR lists the points at which B crosses the real q axis; the largest of these 
is Qc for the given family Gg- The notation "none" in this column indicates that B does not cross the real q 
axis. Column labelled "SN" refers to whether B has support for negative Re{q), indicated as yes (y) or no 
(n). 



Gs 


Ly 


BCy 


BCx 


Nx 


eqs. 


BCR 


SN 


tri 


2 


F 


(T)P 


4 


{2(1),1(2)} 


3, 2, 


n 


**tri 


3 


F 


(T)P 


10 


{3(1),2(2),1(3)} 


3, 2, 


n 


**tri 


4 


F 


P 


26 


{1(1),2(4),1(8),1(9)} 


3.23, 3, 2, 


y 


tri 


3 


P 


F 


1 


{1(1)} 


none 




tri 


4 


P 


F 


2 


{1(2)} 


4, 3.48 


n 


**tri 


5 


P 


F 


2 


(1(2)} 


3.28, 3.21 


n 


**tri 


6 


P 


F 


5 


(1(5)} 


3.25 


n 


tri 


oo 


P 


F 






4, 3.82, 


y 


**tri 


3 


P 


P 


11 


{5(1),3(2)} 


3.72, 2, 


n 


**tri 


3 


P 


TP 


5 


{5(1)} 


3.72, 2, 


n 


tri 


2 


F 


F 


1 


{1(1)} 


none 




tri 


3 


F 


F 


2 


{1(2)} 


2.57 


n 


tri 


4 


F 


F 


4 


{1(4)} 


none 


n 


**tri 


5 


F 


F 


9 


{1(9)} 


3 


n 
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of P and W and studies of B were also presented for wider strips, including (e) cyclic, Ly = 4, 
(f) {PBCy, FBCx) = cylindrical, Ly = 5, 6, and an illustrative {FBCy, FBC^) = open case 
with Ly = 5. A comparative analysis of these results was included. An interesting result of 
our calculations of W on infinite-length strips with different widths and transverse boundary 
conditions is the observation that for the cases studied, B passes through g = 2 (as well as 
g = 0) for strips with periodic or twisted periodic longitudinal boundary conditions but does 
not for strips with free longitudinal boundary conditions. Hence, in particular, if one uses 
periodic or twisted periodic longitudinal boundary conditions, the locus B nicely signals the 
existence of the zero-temperature critical point of the Ising antiferromagnet on these infinite- 
length, finite-width strips of the triangular lattice. Numerical values of W were given for 
infinite-length strips of various widths and were shown to approach values for the 2D lattice 
rapidly. Some conjectures for the behavior of B for arbitrarily wide strips, and for a region in 
the q plane free of chromatic zeros, were stated. These exact calculations of the T = Potts 
antiferromagnet partition function and ground state degeneracy on strips of the triangular 
lattice give valuable analytic knowledge of properties of Potts antiferromagnets. 



Note added: The original version of this paper, submitted in early Oct. 1999, contained 
calculations on Ly = 3 strips of the triangular lattice. In response to a request by a referee 
to perform calculations for wider strips, we have added the results on the Ly = 4 strip 
with {FBCy, PBCx:), the Ly = 5,6 strips with {PBCy, FBC^), and the Ly = 5 strip with 
iFBCy,FBC,). 



Acknowledgment: The research of R. S. was supported in part by the U. S. NSF grant 
PHY-97-22101. 

9 Appendix 

9.1 Generating Function for Ly = 3 Mobius Strip of the Triangular 
Lattice 

As noted in the text, it is convenient to leave the coefficients CtsMbj, j = 6, 7, 8 in the general 
form ( 2.37|) . For the evaluation of these coefficients, we list here the generating function for 
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this strip. We have dj^ = 8 and dx> = 10, and 



10 



V{trt{Ly = 3),FBCy,TPBC^,q,x) = - Xt3j{q)x) 



(9.1.1) 



i=i 



where the Aoj were given in eq. ( |2.15| ). Since several of the Afsj's are algebraic, it is useful 
to display the denominator in an explicitly polynomial form: 



V{tri{Ly = 3), FBCy, TPBC^, q,x) = {l+x) I - {q - 2)x 



1 - (2g - 7)x + {q- 2)(g - 3)x^ 1 + (g - 2)(g - 3)x 



Ffiq'iFfiq2 



where 



i=6,7,8 



(9.1.2) 
(9.1.3) 



with 6^3 J, j = 1, 2, 3 given in eqs. ( |2.18| )-( p^20| ), and 



FtM2 = n (1 - ^*'^-^) = 1 - - + 18g - 17)x + (g - 2)3(g - 3)^^ . (9.1.4) 



i=9,10 



For the numerator, extracting a common factor via the definition of the reduced coefficients 



A 



t3Mb,j^ 



we have 



A 



t3Mb,j 



g(g-l)(g-2)(g-3)A 



t3Mb,j ; 



(9.1.5) 



(9.1.6) 
(9.1.7) 
(9.1.8) 
(9.1.9) 



A3Mb,o = - 6g + 10 

At3Mb,i = {q- 3)(g' - 12g2 + 45g - 55) 

AnMb,2 = q^ - 20g^ + 170g^ - 779g^ + 2016g^ - 2779g + 1588 

AaMb,3 = -{q - 3)(5g^ - 82g^ + 574g^ - 2185g^ + 4745g^ - 5536g + 2691) 

AtsMbA = {q-2){q- 3)(9g^ - 141g^ + 930g^ - 3303^^ + 6651^^ - 7176g + 3224) (9.1.10) 

At3Mb,5 = -{q - 2)'(g - 3)3(6g" - 55g=^ + 186?^ - 277g + 152) (9.1.11) 

AsAfb.e = -{q - 2f{q - 3)(g^ - 12g3 + 59^' - 138g + 125) (9.1.12) 

At3Mbj = {q- 2fiq - 3f{3q^ - Uq + 17) (9.1.13) 

At,Mb,s = -{q-2fiq-3f. (9.1.14) 

It should also be noted that there is a significant difference between the strips of the 
triangular lattice studied here and the analogous strips of the square lattice [|lT], |l^, In 
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general, a strip of the square lattice of width Ly and length with any of the boundary 
conditions {FBCy,PBC^) (cychc), {FBCy,TPBC^) (M5bius), {PBCy,PBC^) (torus), or 
{PBCy, TPBCx) (Klein bottle) is invariant under a translation by one edge or vertex in the 
longitudinal direction. However, in the case of the analogous strip of the triangular lattice, 
this is only true of the cases with cyclic and torus boundary conditions; the strips with 
Mobius and Klein bottle boundary conditions have a "seam" along with the orientation of 
the triangles reverses. This is discussed further in the appendix. Thus, if one proceeds in 
a longitudinal direction along the triangular-lattice Mobius strip, starting in a manner such 
that the triangles are formed by edges connecting the upper left and lower right vertices of 
squares (relative to one's direction of motion), then when one crosses this seam, the triangles 
will be formed by edges connecting the upper right and lower left vertices of the squares on 
the strip. Related to this, there are differences in the degrees of various vertices on these 
strips (where the degree of a vertex is defined as the number of edges that connect to this 
vertex). If one avoids the lowest few values of where the strips degenerate, then, in 
general, (i) for the cyclic triangular strip of width Ly and length Lx-, the {Ly — 2)Lx internal 
vertices have degree 6 while the 2Lx vertices on the upper and lower sides have degree 4; 
(ii) for the same strip as in (i) but with Mobius instead of cyclic longitudinal boundary 
conditions, the {Ly — 2)Lx internal vertices have degree 6, the 2{Ly — 1) vertices on the 
upper and lower sides except for those on the seam have degree 4, and, on the seam, the 
external vertices have degrees 5 and 3; (iii) for the same strip as in (i) but with torus or 
Klein bottle boundary conditions, all of the vertices have degree 6. 

9.2 Equations for the Terms in the Chromatic Polynomial for the 
Cyclic Ly = 4: Strip of the Triangular Lattice 

Four of the Xt4.,j, which we label as j = 2, 3, 4, 5, are the same as for the open Ly = 4 strip 
of the triangular lattice. These are solutions to the quartic equation 

+ &t4,l,l^' + bt4,l,2^^ + btA,l,3^ + 6t4,l,4 = (9.2.1) 



where the coefficients were given as &f(4)j = j = 1, ••,4 in eqs. (B.15)-(B.18) of 

The terms Xt^j, j = 6, 7, 8, 9 are solutions to the quartic equation 

+ btA,2,l^^ + &M,2,2^' + &M,2,3e + &M,2,4 = (9.2.2) 

where 

&M,2,i = 4g - 13 (9.2.3) 

fei4,2,2 = 2{3q^ - ISg + 26) (9.2.4) 
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&t4,2,3 = (? - 2)(4g2 - 25? + 38) (9.2.5) 

and 

bu,2A^{q-2f{q-3f . (9.2.6) 
Another set of A/s for 10 < j < 17 are roots of an equation of degree 8, 

8 

e + ^bu,3,ke-' = (9.2.7) 

k=l 

where 

6*4,3,1 = 2(-3g2 + 19g - 31) (9.2.8) 

6*4,3,2 = 15g^ - 186g3 + 867g2 - 1794g + 1385 (9.2.9) 

6*4,3,3 = -20?^ + 366?^ - 278V + 11248?^ - 25425?^ + 30452? - 15080 (9.2.10) 

6*4,3,4 = (? - 2)(15?^ - 334g6 + 3174g^ - 16676?^ + 52294g^ - 97852?^ + 101138g - 44528) 

(9.2.11) 

6*4,3,5 = -(?-2)2(5-3)^(6g^-126g^ + 1076g^-4804g^ + 11861?^-15378? + 8185) (9.2.12) 

6*4,3,6 = iq- 2)\q - Sfiq^ - 25g^ + 216g^ - 868g^ + 1670g - 1246) (9.2.13) 

6*4,3,7 = (g - 2f{q - 3)^(2g3 - 22q^ + 8Qq - 97) (9.2.14) 

6*4,3,8 = (g-2)«(?- 3)^ . (9.2.15) 
A final set of At4j, 18 < j < 26, are solutions to the equation of degree 9 

9 

e + J2buA,ke-''^0 (9-2.16) 



k=l 



where 



6t4,4,i = 2(g - 3)(2?2 - 12g + 21) (9.2.17) 

6*4,4,2 = iq- 3)(6g^ - 90g^ + 558q^ - 1772g2 + 2865g - 1875) (9.2.18) 

6*4,4,3 = - lllg^ + 1380g^ - 10071g^ + 47476g^ - imi2q'^ + 315652g=^ 
-428385?^ + 339300? - 119368 (9.2.19) 

6u,4,4 = {q-2){q- 3)2(g9 - 3Aq^ + 491g^ - 4032?^ + 20961g^ - 71954?^ 

+163654?^ - 238278?^ + 201722? - 75672) (9.2.20) 
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6^4,4,5 = -{q - 2)^(g - 3f{3q^ - 81q^ + 983g^ - 7029g^ + 32609g^ - 101701g^ 
+213036g^ - 288702^^ + 229385g - 81299) (9.2.21) 

buAfi = iQ- 2)^(g - 3)^(3g^ - 70g^ + 728g^ - 4404^^ + 16929g^ - 42286g3 

+66933g2 - 61296g + 24830) (9.2.22) 

6^4 47 = -(g - 2)%q - 3f{q^ - 20q^ + 166g^ - 734^^ + 1833g^ - 2462g + 1393) (9.2.23) 
6^44g = -(g - 2f{q - 3)^(2g^ - 21^^ + 87^^ - 165g + 119) (9.2.24) 
6*4,4,9 = -(g-2)^'(g-3)^ (9.2.25) 

9.3 Generating Function for the Ly = 5,6 Strips of the Triangular 
Lattice with {PBCy, FBC^) 



For the Ly = 5 strip we calculate a generating function of the form ( p.2|) with dj) = 2, dj\j- = 1 
and, in the notation of eqs. (|2.3|) and (2.4), 



bt5PF,2 = {q- 3)(3g^ - 66g^ + 619g^ - 3205g^ + 9877g^ - 18065g^ + 18078g - 7588) (9.3.2) 



^5PF,o = q{q - l)(g - 2)(g - 3)(g^ - 14g^ + 85g^ - 290g^ + 599^^ - 723g + 398) (9.3.3) 

At^PF,! = -q{q - l)(g - 2)(g - 3){q^ - 2q + 2)(3g^ - 66g^ + 619g^ - 3205g^ 

+9877g^ - 18065g2 + 18078g - 7588) . (9.3.4) 

For the Ly = 6 strip we calculate a generating function of the form ( p.2|) with dx> = 5, 
d^f = 4 and 

btQpp^^ = -q^ + 18g^ - 145g^ + 680g^ - 1980?^ + 3379g - 2586 (9.3.5) 
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h6PF,2 = 4g^° - 128q^ + 1868g^ - 16352g^ + 94977g^ - 382031g^ + 1076317g^ 
-209389V + 2686606g^ - 2047842? + 702080 (9.3.6) 

bt6PF,3 = -2(g - 3)2(gi2 _ 39^11 + QQ^io _ 7493^9 ^ 54599^8 _ 282283g^ + 1068575?^ 
-2982861g^ + 6098756g^ - 8908956?^ + 8820488?^ - 5306146? + 1462992) (9.3.7) 

bt6PF,4 = -2{q -2){q- 3)^{2q^^ - 62q^^ + 882g^ - 7601g^ + 44105g^ - 181018?^ 
+536536g^ - 1149015?^ + 1742334?^ - 1779827?^ + 1099188g - 309188) (9.3.8) 

bt6PF,5 = 4(g - 2)2(g - 3)^(g2 -5q + 5)(2g^ - 36g^ + 277g^ - 1179g^ + 2990g2 
-4505?^ + 3728g - 1310) (9.3.9) 

AePFfi = q{q -!)(?- 2)(g^ - 21^^ + 199g^ - 1121?^ + 4159g^ - 10623g^ 
+18887g^ - 22824g2 + 17177? - 6143) (9.3.10) 

At6PF,i = -qiq - 1)(9 - 2)(4gi^ - 140?^^ + 2266g^^ - 22416g^° + 150973?*^ 

-730186g^ + 2607252g^ - 6958852?^ + 13899608?^ - 20584349g^ 

+22103679?^ - 16461349?^ + 7723994? - 1748140) (9.3.11) 

At6PF,2 = 2q{q -l){q- 2){q - Sfiq^^ - A2q^^ + 816?^=^ - 9734?^^ + 79793g^^ 
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-476549g^° + 2144264g^ - 7409966g* + 19852299g^ - 41297346g^ + 66301130g^ 
-80939629g^ + 73099740g^ - 46448750^^ + 18742947g - 3655548) (9.3.12) 

At6PF,3 = 2qiq - l)(g - 2)(g - 3)^(2g^^ - 72g^^ + 1212^^^ - 12651^^^ ^ 91556^11 
-486599g^° + 1962326g^ - 6116898g^ + 14870220g^ - 28223745g*^ + 41554711g^ 
-46735427g^ + 39026466^^ - 22975509?^ + 8586616g - 1545752) (9.3.13) 

AtGPFA = -Hi - - 2)2(g - 3fiq^ -5q + 5)(g^ - 5g=^ + lOq^ - lOq + 5) 

X (2g^ - 36g^ + 277g^ - 1179g^ + 2990g3 - 4505?^ + 3728g - 1310) (9.3.14) 

9.4 Generating Function for Ly = 5 Open Strip of the Triangular 
Lattice 

For this strip we calculate a generating function of the form ( |2.2| ) with rf© = 9, dj\/ = 8 and, 
in the notation of eqs. (|2.3|) and ( p.4|) , 

bt5FF,i = -{q - 3)(g' - lOg^ + 46q^ - 113g + 120) (9.4.1) 

bt^pp2 = 6q^ - 131g^ + 1280g^ - 7328g^ + 26930g^ - 65081?^ + 100888?^ - 91462g + 36965 

(9.4.2) 

5^5^^3 = _(g - 2)(g - 3)(15g^ - 378g^ + 4289g^ - 28788g^ + 126096g^ - 374139g^ 
+752541^3 - 989867g^ + 77261 Ig - 272483) (9.4.3) 

bt5FF,4 = (g - 2)(g - 3)2(20g" - 607g^° + 8429g^ - 70702g^ + 398115g^ 
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-1580547g^ + 4515585?^ - 9285872?^ + 13471537g=^ - 13131321g^ 

+7738560g - 2087938) (9.4.4) 

bt5FF,5 = -(q - 2)^(g - 3)3(15?^^ - 502?^^ + 7729g^° - 72397?^^ 
+459566g^ - 2083176g^ + 6915864g^ - 16947196g^ + 30430188g^ 
-39053679g^ + 34008163g2 - 18041392g + 4408580) (9.4.5) 

bt5FF,e = (? - 2)^(? - 3)^(6g^2 - 203g^^ + 3149?^° - 29630?^^ + 188440?^ - 853749?^ 

+2826657g^ - 6893466g^ + 12293272g^ - 15636918?^ + 13466958g^ 

-7049578g + 1695556) (9.4.6) 

bt5FF,7 = -(g - 2)^(g - 3)^(g^^ - 34g^° + 516g^ - 4647g^ + 27718?^ - 115308?^ 
+342008g^ - 724072g^ + 1072964g3 - 1060043?^ + 628196g - 169014) (9.4.7) 

bt5FF,8 = -{q - 2)^(g - 3)^(g^ - 23q^ + 241g^ - 1505?^ + 6145g^ - 16929?^ 
+31319g^ - 37359?^ + 25972g - 7987) (9.4.8) 

W,9 = (g - 2y\q - 3y\q' - 8q^ + 21q - 17) (9.4.9) 
We use the definition At^ppj — q{q — l){q — '2,)At^Fpj and have 

At5FF,o = {q- 2)' (9.4.10) 
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At^FF,i = -{q - 2f{6q^ - 113g^ + 926g^ - 4304g^ + 12381g^ - 22504g3 
+25133?^ - 15663? + 4121) (9.4.11) 



At5FF,2 = {q- 2f{q - 3)(15g^° - 363g^ + 3936g^ - 25141g^ + 104572?^ 
-295342g^ + 572184g^ - 748554?^ + 630336?^ - 306798? + 65006) (9.4.12) 

At5FF,3 = -{q - 2f{q - 3)2(20?^^ _ 537^11 + 7867^10 _ 63587?^ + 344801?^ 
-1319670?^ + 3650087?^ - 7338694?^ + 10623070?^ - 10770336?^ + 7237158?^ 
-2882235? + 511741) (9.4.13) 

At5FF,4 = iq- 2)^(? - 3)'^(15?^^ - 487?^2 + 7266?^^ - 65902?^° + 405008?^ 
-1779077?^ + 5739974?^ - 13753341?^ + 24433426?^ - 31724357?^ + 29191727?^ 
-17973781?^ + 6613085? - 1093258) (9.4.14) 

At5FF,5 = -{q - 2)^(? - 3)^(6?^^ - 215?^^ + 3554?^^ - 35907?^^ + 247594?^° 
-1231849?^ + 4556733?^ - 12718548?^ + 26883798?^ - 42760478?^ + 50287310?^ 
-42307379?^ + 24007733?^ - 8200007? + 1267630) (9.4.15) 

At5FF,6 = {q- 2)^(? - 3)^(?^2 - 33?^^ + 485?^° - 4225?^^ + 24378?^ - 98266?^ 
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+283651g^ - 589907g^ + 875195g^ - 900611^^ + 607987g^ - 240848g + 42203) 

(9.4.16) 

Atr^FFj = {q- 2)^°(g - 3)^(g^ - 18g^ + 145g^ - 674g5 + 1941g^ - 3474g^ 
+3701g2 - 2120g + 499) (9.4.17) 

At5FF,8 = -{q - lf{q - 2)"(g - 3)"(g3 - 8^^ + 2lq - 17) . (9.4.18) 

9.5 Ly = 2 Cyclic and Mobius Strips of the Triangular Lattice with 
Odd Nt 

In this part of the appendix we shall report some new results for Ly = 2 cyclic and Mobius 
strips of the triangular lattice with an odd number Nt of triangles and compare these with 
the case of even A^^^ = L^Ly. Some illustrative strip graphs are shown in Fig. ||. 




(c) (d) 

Figure 8: Illustrative strip graphs of the triangular lattice with width Ly — 2: (a) cyclic, Lx — 4, even 
Nt = 2Lx] (b) Mobius, = 4, even Nt = 2L^; (c) cychc, odd N = 7; (d) Mobius, odd Nt = 7. 

For the Ly = 2 cyclic strip of the triangular lattice with even Nt = 2Lx we calculated 
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|n[ (see also 0, P) 

Pitrt{Ly = 2, eye), N, = 2m, q) = {q^ - 3g + 1) + [(g - 2)Y + (?-!) (^^2,3)™ + {^2,4)' 



where 



1 r 



(9.5.1) 
(9.5.2) 



Ai2,(3,4) = ^ [5 - 2g ± y/9-4q 

For the corresponding Ly = 2 Mobius strip of the triangular lattice with even Nt = 2L^, we 
calculated 

(At2,3)™ — (At2,4)' 

P{trz{Ly = 2, Mb.), Nt = 2m,q) = -1 + [{q - 2)Y -(?-!)(?- 3)^ ^ ^ 

'^t2,3 ~ ^t2,4 

(9.5.3) 

Note that the Mobius strip has a seam. The function W and the boundary B were given in 



n[; B separates the q plane into three regions and crosses the real axis at g = 0, 2, and at 



Qc = 3. 

We proceed to our new results. For the Ly = 2 cyclic strip of the triangular lattice 
containing an odd number A^^ = 2m + 1 of triangles, we calculate 



P{trt{Ly = 2, eye), Nt = 2m + 1, q) = -(g - 3g + 1) + (g - 2)[(g - 2)T+ 



1 



(g-l)(g-3) 



(Ai2,3)'" + (At2,4)'" + 



(A 



t2,3) 



(At2,4)' 



(9.5.4) 



At2,3 ~" At2,4 

This strip has a seam. 

For the Ly = 2 Mobius strip of the triangular lattice with an odd number Nt = 2m + 1 
of triangles (which does not have a seam) we obtain 

P{tm{Ly = 2, Mb.), Nt = 2m + l,q) = l + {q-2) [(g - 2)Y+ 



(Ai2,3)"^ + (Ai2,4)'" + (9 - 4g) 



(At2,3)'" "~ (At2,4)^ 



(9.5.5) 



At2,3 ~ At2,4 

Note that the sum of the coefficients C = for the chromatic polynomials for both of the 
Ly = 2 strips with odd A^^, eqs. ( 9.5.4 ) and ( |9.5.5 ). For comparison, for the Ly = 2 even-Nt 
strips, C = g(g — 1) for the cychc case, eq. ( p.5.1| ) and (7 = for the Mobius case, eq. ( p.5.3| ). 

Equivalently, one may write these results in terms of generating functions, with the 
definition analogous to (pTll): 

00 

r{tri{Ly = 2, BC^), Nt = even, g, x) = ^ P{tri{Ly = 2, BC^), Nt = 2m, g)x""^ . (9.5.6) 

m=2 



47 



and 



T{tri{Ly = 2, BC^), Nt = odd, g, x) = ^ P{tri{Ly = 2, BC^), Nt = 2m + 1, q) 



m-2 



X 



(9.5.7) 



m=2 



where BC^ = PBC^ or TPBC^- The generating functions for all of these strips have the 
same denominator, 



V{tri{Ly = 2, eye.)) = JJ(1 - Xt2j 

1 - (g - 2)^x] [l - (5 - 2g)a; + (g - 2)^x^ 



1 -X 



(9.5.8) 



The numerators are easily worked out from the results that we have given for the chromatic 
polynomials and the denominator V; for example, 



T{tri{Ly = 2, eye.), Nt = even, q, x) 



T(tri{Ly = 2, eye.), Nt = odd, q, x) 



q{q - l)(g - 2) g - 3 + (g - 2)x - (g - 2)V 



V{trt{Ly = 2,eye.)) 
g(g-l)(g-2)(g-3)[g-3 + (g-2)2x] 



(9.5.9) 
(9.5.10) 



V{trt{Ly = 2,eye.)) 
and so forth for the Mobius strips. 

We remark on some general features of these results. The chromatic polynomials for both 
even and odd Nt and both cyclic and Mobius strips have the same four Aj's, and hence the 



same W functions and boundary B (given in ITT]]). For all cases, gc = 3. These properties are 



in agreement with the general discussion in |TT|, |T^, ^ on the effects of boundary conditions 
on W and B. 
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